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Two remarkable features emerge from the exact Wilsonian procedure for integrating out the high- 
energy scale in the Hubbard model. At low energies, the number of excitations that couple minimally 
to the electromagnetic gauge is less than the conserved charge, thereby implying a breakdown of 
Fermi liquid theory. In addition, two charge e excitations emerge in the lower band, the standard 
projected electron and a composite entity (comprised of a hole and a charge 2e bosonic field) which 
give rise to poles and zeros of the single-particle Green function, respectively. The poles generate 
spectral weight along an arc centered at (7r/2, 7r/2) while the zeros kill the spectral intensity on the 
back-side of the arc. The result is the Fermi arc structure intrinsic to cuprate phenomenology. The 
presence of composite excitations also produces a broad incoherent pseudogap feature at the (tt, 0) 
region of the Brillouin zone, thereby providing a mechanism for the nodal/anti-nodal dichotomy 
seen in the cuprates. 



I. INTRODUCTION 

As revealed by extensive angle-resolved photoemission 
(ARPES) studies^ lightly doped copper-oxide super- 
conductors (cuprates) in the pseudogap regime possess a 
band of excitations that only cross the chemical potential 
once. Such a single crossing generates a set of coherent 
or pole-like excitations that ultimately form a truncated 
Fermi region, termed a Fermi arc, as opposed to the tradi- 
tional Fermi surface generated by a double crossing. The 
coherent excitations, centered around (0,0) (7r,7r), 
traverse the zone diagonal and terminate in the vicinity 
of (tt, 0) or (0,7r), thereby giving rise to a nodal/anti- 
nodal dichotomy^^^^, the former being ungapped while 
the latter is gapped. While some ARPES experiments^^ 
performed on Bi2Sr2-ajLaa^Cu06+5 (La-Bi2201) revealed 
closed hole pockets, and hence consistency with the re- 
sults from quantum oscillation experiments in high mag- 
netic field^— , this interpretation has been called into 
question^. King, et ai^ observed that the closed pock- 
ets seen earlier— are entirely structural in origin as they 
originate from overlapping superstructure replicas of the 
main and shadow bands. Consequently, the preponder- 
ance of evidence from ARPES is that the coherent exci- 
tations form a disconnected region in momentum space 
consistent with a single crossing of the chemical poten- 
tial. 

Theoretically, two questions arise. 1) What suppresses 
the spectral weight on the back-side of the arc? 2) What 
is the origin of the incoherent excitations or gap at the 
zone boundary? A natural candidate to explain the for- 
mer is that two kinds of excitations populate a doped 
Mott insulator, one of which has no overlap with the 
electron. Such an excitation will appear in the single- 
electron Green function as a zero rather than a pole and 
hence will carry no spectral weight. In this vein, some 
have proposed neutral composite excitations^^ to explain 
the origin of Fermi arcs. Alternatively, Fermi arcs have 
been seen in numerics^Sri^^ on the 2D Hubbard model 
and have been modeled phenomenologically^^ (hereafter 



YRZ). However, a key assumption of the phenomenolog- 
ical account is that the zero-line is fixed at the diamond- 
shaped Fermi surface of the non-interacting system. That 
the diamond-shaped Fermi surface of the non-interacting 
system constitutes the zero-line of the single-particle 
Green function is a rigorous mathematical statementi24 
only if the underlying Hamiltonian is particle-hole sym- 
metric. In fact, from the precise condition2i~i2£ for the 
vanishing of the real part of the Green function, main- 
taining that the zero-line is doping independent requires 
unphysical assumptions regarding the spectral function. 
Certainly such a conservation of the zero-line is not 
borne out by numerics on the Hubbard model2ii2^ nor 
by analytical arguments^^. In addition, models involving 
Cooper pairs, fluctuating or otherwise, have been con- 
structed to either yield arcs^ or hole pocket o^^i^^ . Our 
discussion here, however, will focus entirely on arcs as 
this seems to what the ARPES experiments are about. 

While the physical origin of arcs might not be clear, 
the mathematics is. Any Green function of the form 
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will do. The similarity with the BCS Green function is 
only perfunctory as there is no anomalous component. 
Aside from having poles, Eq. (jl.ip has zeros whenever 
= assuming of course the chemical potential cor- 
responds to a; = 0. While the two dispersing electronic 
bands, ek and ^k in Eq. (|l.ip are not hard to come by, 
the parameter A(k) is. It requires some sort of order, 
fluctuating or otherwise, or a new bosonic degrees of 
freedom. While YRZ proposed Eq. (jl.ip phenomenolog- 
ically, their intuition was based on weak-coupling RPA 
diagramatics^^ on 2-leg ladder systems. In the context 
of an algebraic charge liquid (ACL)^^, Qi and Sachdev^^- 
also obtained Eq. (jl.ip . Thus far, the only system shown 
to have the properties of the ACL is one with radically 
different parameters than the cuprates^. 

Our point here is that A(k) arises fundamentally from 
a new degree of freedom associated with dynamical spec- 
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tral weight transfer. Ideahy, it would be advantageous to 
derive Eq. (|l.ip directly from the strong-coupling param- 
eter space of the basic model for a doped Mott insulator, 
for example the Hubbard model. Such a derivation has 
not presented in the literature. Hence, it is this problem 
that we address. Since zero a^^i^^ arise from a cancellation 
of the spectral weight in the upper and lower Hubbard 
bands (hereafter UHB and LHB), an accurate description 
of the upper band is required in a derivation of Eq. (jl.ip . 
Hence, an attempt to obtain Eq. (|l.ip from a model that 
projects out the UHB, for example the t — J model, is 
initially a non-starter as this model does not have zeros 
of the type required for Eq. (jl.ip . Nonetheless, the infor- 
mation regarding the UHB should be correctly encoded 
into a theory of the lower band if the UHB is integrated 
out (rather than projected out) exactly ala Wilson^^. In 
this paper, we show how the method we have recently 
developed^— for carrying out the Wilsonian program 
for the Hubbard model can be used to derive Eq. (jl.ip . 
We show expHcitly that two types of excitations emerge, 
projected electrons (yielding poles in the propagator) and 
a new bound state that gives rise to zeros. The bound 
state is not made out of the elemental excitations and 
hence is orthogonal to an electron (hence the zero) . It 
represents a charge e excitation that originates from the 
non-rigidity of the Hubbard bands, in other words, the 
well documented dynamical spectral weight transfer— —, 
the key fingerprint of the breakdown of the band concept 
in Mott systems. Since the mathematics of Fermi arcs 
requires two kinds of excitations, one with poles and the 
other with zeros, we refer to a physical model that con- 
tains both, such as the one presented here, as the stan- 
dard model. More exotic models relying on some type 
order would fall outside this framework. 



II. CHARGE 2e BOSON THEORY 

A. Preliminaries 

Underlying our toy Green function is a two-pole struc- 
ture of the form, 

^ / , x r. f cos^ sin^ \ , , 
Gtoy(a;, k) = Z + . (2.1 

\UJ — UJj^ UJ — UJ- J 

Here cos^ 6 = (cj+ — ^(k))/ (cj_ — cj+) and uj± = ^ (e(k) + 
£:(k)) ± v'l^W -^(k))2 +4|A(k)|2. Zeros arise from the 
interference between the poles at uj- and ujj^. Any model 
that admits zeros must have at least this two-pole struc- 
ture. There are two limits of the Hubbard model in which 
the zero surface can be calculated exactly. In the atomic 
limit, the zero surface of the exact single-particle Green 
function 

G«(c.) = + (2.2) 

is independent of momentum given hy uo = (i and x = 0. 
When the hopping is non-zero, t 7^ 0, the only limit in 



which the zero surface can be calculated exactly is at 
half-filling and particle- hole symmetry. In this limit, the 
zero surface^ for a nearest-neighbour band structure is 
the magnetic Brillouin zone. Since Fermi arcs are absent 
from both the atomic limit and the half-filled system with 
hopping, it follows necessarily that Fermi arcs (if they are 
present at all in the Hubbard model) arise entirely from 
the dynamical part of the spectral weight. 

Dynamical spectral weight transfer represents a con- 
crete example of more being totally different. As is evi- 
dent from Eq. ()2.2p . the weight of the lower band in the 
atomic limit is 1 + x. This spectral weight has a natu- 
ral interpretation in terms of electron states. There are 
2x electron addition and 1 — x electron removal states. 
Hence, in the atomic limit, there is a one-to-one corre- 
spondence with the spectral intensity and the number of 
electron states in the lower band. When the hopping is 
turned on, the spectral intensity increases in the lower 
band to 1 + x + a, where 

"=|^E</U.) + 0((W) (2.3) 

and the /icr's are a rotation of the original fermionic oper- 
ators in the Hubbard model such that it is block diagonal. 
The energy of each block is n/7, n the number of dou- 
ble occupancies in each block, a is necessarily positive 
because any hopping process that creates doubly occu- 
pancy decreases the available spectral weight in the upper 
band. Consequently, counting electrons, fractionalized or 
otherwise, cannot exhaust the total number of degrees of 
freedom in the lower band. A new degree of freedom 
must be present which is distinct and hence orthogonal 
to electron quasiparticles. This degree of freedom will 
appear as a zero in the spectral function in the lower 
band. It is precisely the nature of the states that arise 
from the mixing with the upper band that we elucidate 
here. 

Several approaches suggest some kind of composite ex- 
citation mediates Fermi arcs. Consider the SU (2) gauge 
theory of the t-J model proposed by Wen and Lee^ 
in which the elemental fields are the appropriate linear 
combinations of two charge e bosons and two spinons. 
A mean-field calculation^ of the coherent spinon-boson 
Green function reveals that the spectral weight of the oc- 
cupied part of the spectrum exceeds 1 — x, acquiring a 
value of 1 + x/2 instead. Wen and Lee^ alleviated this 
problem by introducing an interaction which recombined 
the bosons and fermions back into the elemental fields. 
The effect of this interaction with strength U was to en- 
hance the spectral weight in the unoccupied part of the 
spectrum. The correct spectral weights were obtained 
simply by adjusting the magnitude of U. The unnoc- 
cupied part of the spectrum^^ shows up as a small hole 
pocket centered roughly at (7r/2,7r/2). The intensity on 
the back-side of the pocket is greatly suppressed, thereby 
leading to a structure not too distinct from a Fermi 
arc. Subsequent work on the U{1) formulation of the 
t-J model with a phenomenological spinon-holon binding 
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term reached a similar conclusion^. However, the most 
extensive calculation in the gauge theoretic formulations 
of the Hubbard model reached a rather different conclu- 
sion. Working directly with the Hubbard model, Imada 
and colleagues^^ used a slave-particle construction with 
the gauge fluctuations treated at the RPA level and con- 
cluded it is actually dynamical spectral weight transfer 
that leads to a fermi arc structure and a Green function 
of the YRZ form. Aside from suffering from the lack of 
a systematic way of treating the gauge fluctuations, this 
formulation generates Fermi arcs from a neutral compos- 
ite excitation^. Since neutral entities cannot couple to 
the current, it is unclear how such entities can influence 
the spectral function. What we demonstrate here with- 
out resorting to a gauge theory is that dynamical spectral 
weight transfer mediates Fermi arcs. 



B. Exact Results: The conserved charge and the 
low energy mode 

Dynamical spectral weight transfer has two profound 
consequences. First, we show exactly that the conserved 
charge {1—x) is not exhausted by counting the degrees of 
freedom minimally coupled to the electromagnetic gauge. 

I 



The remainder are carried by an incoherent background. 
Second, the dynamically transferred degrees of freedom 
give rise to zeros in the lower band. In fact, the physics 
we find here is analogous to the finite U charge pro- 
cesses that contribute to the spectral function measured 
by ARPES in heavy-fermions^. 

We first start with the procedure to integrate out the 
UHB. Previously, we demonstrated^^^^ that a theory of 
the lower Hubbard band (hole doping) can be obtained 
by introducing a new fermionic field Di^^ which rep- 
resents the excitations in the upper band. This field has 
mass U and hence should be integrated out to construct 
the exact low-energy theory. As in the derivation of all 
collective phenomena which dates back to the classic pa- 
per of Bohm and Pines^^, a constraint must be intro- 
duced so that when solved, the model in the extended 
space is equivalent to the starting UV-complete theory. 
The new fermionic degree of freedom enters the action 
in a quadratic fashion and hence the standard fermionic 
path integral techniques can be applied to integrate the 
high-energy scale of the UHB. Since the lower and upper 
bands are not rigid in the sense that the spectral weights 
of the two bands are coupled, integrating out the upper 
band will lead to new degrees of freedom in the lower 
band. In Euclidean signature, the Hubbard action in the 
extended space is 



- [^^(1 - - nja)cl^Cja + D*c*^Ci^Dj + {D^Oci^V^Cj^ + c.c.) + se'^i^*{Dj - Ocj^cji) + c.c. I, 

(2.4) 



where the matrix gij selects the relevant neighbors, = 
±l(cr I), the constraint is given by S{Di — Oci^Cii)^ 
^ is a Grassmann, 5 is a constant with units of energy 
so that cpi is dimensionless and q^t is an electron anni- 
hilation operator for site i with spin a. Because the S- 
function constraint appears exponentiated in the action, 
an auxiliary field with charge 2e, cpi enters the action 
as a Lagrange multiplier. As a consequence, the field 
(f is not made out of the elemental excitations (thereby 
distinguishing it from other charge 2e scenarios involv- 
ing pairs of electrons) but rather arises because the UHB 
and LHB are not rigid bands. In the action, the first 
two terms represent the dynamics in the lower and up- 
per Hubbard bands, respectively, the third term the mass 
of the D field, the fourth term the hopping in the lower 
band with matrix element the next two the dynamical 
mixing between the upper and lower bands and the last 
term the constraint. The constant s has units of energy 
and is 0(t)^. It is straightforward to check that solving 



the constraint by integrating out the auxiliary field, ipi^ 
followed by an integration over Di exactly reduces 
to the action for the standard Hubbard model. This is 
the UV limit of our theory. The advantage of the refor- 
mulation above is that it cleanly associates the physics 
of the upper band with a fermionic field Di which enters 
the action in a quadratic fashion. To obtain the IR limit, 
one simply has to perform the Gaussian integration over 
the massive field, D^. The result is the low-energy or IR 
action, S]^ = J drC]^, with the associated Lagrangian, 

^h^ = (1 - nia)Ci^drCi^ - tgij(l - ni^)cl^Cj^{l - rij^) 

- {sifi - thiY {M-^).. {sLpj - thj) 

- (s(/:?*Q^C4 + C.C.) - -^trlnA^, (2.5) 

where a matrix element of M is given by Mij = {dr + 
U)8ij - tgijC^j^Cia- and bi = QijCjcjVcjCi^. Hereafter, 
repeated indices are implicitly summed unless otherwise 
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stated. It is important to note that no approximations 
have been made as of yet. 

In both actions, S}^^ and Sj^^ global U{1) symmetry 
guarantees the existence of a conserved charge, which 
turns out to be 

gr = (l-n.^)c*,c,, + 2D*A, (2.6) 

Q/i^ — ~ ^id-)CicF^io- 

+ 2 {s^i - tbif (A^-i) {M-')^. {s^j - tbj) . 

(2.7) 

It is a natural consequence that the conserved charge 
{Q]l^) is consistent with the number of electrons in the 
original Hubbard model since the operator DjDi counts 
the number of doubly occupied sites. From the Hellman- 
Feynmann theorem, it is straightforward to check how 
the number of double occupancies, ^^occ^^' related 
to A. Since ri^^JJ^ = p-^d\n zj^"^^^^ /dU with Zh = 
J ]e~*^^, one can easily observe the second terms 
in Eq. (|2.6p and (|2.7p are identical to n^^^ and ^docc 
respectively. As a result, the conserved charge {Q]^) is 
identified with the number of electrons, 1 — x, with x 
the number of holes. This is one of the indications that 
the low-energy action, Sj^, retains the structure of the 
Hubbard model even after the integration of the massive 
modes. 

Another advantage of the low energy theory is that the 
non-Fermi-liquid nature of the low-energy excitations is 
immediately manifest. To illustrate this, one can add a 
minimally coupled source term^^ 

= Jl [^^(1 - nis)ci„ + ec*aVaDi] + c.c. (2.8) 

SO that when the constraint is solved, the bare electron 
operator is generated^^. What we would like to know is 
what is the transformed fermion at low energies. Inte- 
grating out the Di fields, results in a source contribution 
to the IR Lagrangian, 

C'f = Jl^P,, + c.c. - JIC*-, {M-^)^^ Cj-aJja (2.9) 

with a new collective field, ^jjia- given by 

= (1 - n,,)c*, + tb* {M-').. V.Cis 

-s^*(M-i)..Kci^. (2.10) 

Note V^icr is derived not contrived. It is Eq. (9) in the 
Ref. 35. We obtained it by integrating the UV— com- 
plete Lagrangian in the presence of the source term that 
generates the correct UV current with respect to the mas- 
sive field Di. tjjia- is the propagating degree of freedom 
in the IR. It contains not only an electron-like quasipar- 
ticle affected by nearby spin fiuct nations, but also a hole 
(with the opposite spin) that is dressed with a doubly- 
charged bosonic mode. Note, we cannot give tl^ia- a simple 
interpretation in terms of bosons or fermions. At best, 
ipia- corresponds to the physical field that is minimally 



coupled to an external gauge field. That is, these are the 
excitations that couple to light. Hence, it is the field that 
is probed by an ARPES experiment, for example. While 
i/jia- was derived earlier, what we did not show explicitly 
is that it does not stand in a one-to-one correspondence 
with the bare electrons. This can proven exactly by fo- 
cusing on the positive-definite quantity, 

^ia^ia = (1 - ^m)c*aCi(7 

+ {tb - s^Yj [M-')^^ c,,c:, [M-^)^^ {tb - s^)^ 

= Qf 

- {tb - s^Tj (2 + c:,c,,) {M-')^^ {tb - s^)^ , 

(2.11) 

which is essentially the conserved charge less the num- 
ber of doubly occupied sites. Since the second term in 
the last line is positive definite, the number of low-energy 
collective modes which are minimally coupled to the elec- 
tromagnetic gauge field is less than {Q]^) = 1 — x. The 
natural resolution of this conumdrum is that the num- 
ber operator only counts those excitations that have a 
particle- like interpretation. That is, the number opera- 
tor only counts the coherent part of the spectrum. All 
of the stuff mediated by mixing with the upper band is 
entirely incoherent and hence while it can contribute to 
the current, it is not enumerated by counting the number 
of particles. The remainder of the charge count is carried 
by the last term in Eq. (|2.9p . 

This discrepancy is not a surprise when one considers 
that the total spectral weight of the lower band exceeds 
1 + Q[S^^ by a dynamical correction, a > 0, that de- 
pends on the hopping integral, t. Since there are only 
1 + X electron states in the lower band, and only charge e 
excitations contribute to the spectral function, there has 
to be some new charge mode to make up the difference. 
What ijjicr lays plain is that there are charge e states that 
contribute to the current that are completely incoherent. 
It is a composite excitation of Lp"^ and a hole ci^. In terms 
of the UV variables, this degree of freedom represents the 
binding of a doublon and a holon. The new composite 
excitation, (/?*A^~^]^Ca- has internal structure and hence 
is orthogonal to the projected electron. Since there is no 
Hilbert space for interpreting (^"^ Cia- in terms of a parti- 
cle is not possible. It is this additional degree of freedom 
that creates the Fermi arc structure-that is, the zeros of 
of the Green function. Hence, hidden in ipia- is an inco- 
herent contribution to the single-particle Green function. 
What this discussion makes clear is that should not be 
considered to be an independent degree of freedom but 
rather one that is strongly coupled to the fermions. 

What we have shown thus far is that there is a dy- 
namical contribution to the charge degrees of freedom 
that are coupled to the source term that generates the 
current. Such entities are the physical degrees of freedom 
that create holes in the lower band. Consequently, when 
one such excitation is removed from the lower band, the 
change in the spectral weight should also depend on t. 
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Hence, the doping level should receive a dynamical con- 
tribution. To this end, we defined^^ x' = x + a and 
hence the weight in the UHB is 1 — and the occupied 
and empty parts of the lower band are 1 — x' and 2x' ^ 
respectively. 

C. Green function and the approximations 

Thus far, all of our statements are exact. Our calcula- 
tion of the Green function is not, however. To lend cre- 
dence to our treatment, we state our assumptions clearly 
and up front. The complexity arises in treating the ip 
degree of freedom. Our treatment is in the spirit of the 
results obtained in the previous section, namely that ip 
leads to the creation of a new charge e excitation that 
is orthogonal to a projected electron on account of its 
internal structure. 

Having determined the generating functional, 
^' h[{Ji(ji Jicr}]i we proceed to calculate the Green 
function. In the functional formalism, it is given by 

= - {TrMr)^]{0)) + ((5(r)4 {M-')^^ c^,) , (2.12) 

where % represents time ordering and (• • •) stands for the 
average over all possible paths. Since the second term is 
independent of time, this term contributes to the inco- 
herent part of the Green function. To understand the 
first term which contains both coherent as well as inco- 
herent responses, it is sufficient to focus on the correlator 
between the ?/^icr's. Since i/jia contains a composite excita- 
tion which has a prefactor of t/U, {Trivia {^)i^jcr{^)) can, 
in principle, be expanded in power oft/U. The presence 
of the projection operator, (1 — n^a^), however, does not 
necessarily guarantee that it gives a dominant contribu- 
tion compared to the composite entities. 

For the purpose of numerical evaluation, we make an 
approximation to the projection operators, following the 
idea developed by Zhang, et alA^ However, the crucial 
difference is that we make the substitution the bare hole 
concentration by the effective hole doping level (x ^ x'), 
since the physical entities coupled to the external gauge 
field are not the bare electrons but are rather dynamically 
generated. Hence our first approximation is 

(A-1) (1 - ni^)c*^c^-^(l - rij^) gtC^^Cj^, 

where gt = 2x^/(1 -\-x^). Interestingly, in the strong cou- 
pling limit, {U/t 1), di mean- field approach to Kotliar- 
Ruckenstein's slave boson construction^^ led to the same 
renormalization factor for the charged fermion but with 
x^ replaced with x. Likewise, 

(A-2) (1 - ni^)cl^{dr H ^ dpclai^r H )Cia- 



where gp = {1 — x') / {1 — x). The multiplicative factors 
are chosen here for internal consistency with the two as- 
sumptions. 

Since the action has all relevant degrees of freedom 
for the low-energy sector, including the spin singlet fluc- 
tuations {bi) and mixing between the separate Hubbard 
bands ((/?), it is reasonable to expand the action in powers 
oit/U. To leading order, the matrix elements {M~^)ij is 
U~^5ij. From the fact that the collective boson (fi only 
has dynamics through its coupling to the fermions, we 
assume the dynamics of the boson to be frozen. Opera- 
tionally this assumption breaks down at 0{t/U)'^ where 
the explicit dynamics of (f appears as can be seen from 

2 

an expansion of the A4 matrix, — U + • • • 

In fact, this even at 0{t/U)'^, the propagator for cp lays 
plain that it has a pole only in the high-energy sector. 
This justifies the assumption that 

(A-3) Bosonic field, has no dynamics in the LHB. 

In other words, it alone is highly massive and is not likely 
to propagate in the low-energy sector. Finally, although 
local spin ordering might non- negligible, we will assume it 
to be at most ancillary to the strong interaction physics 
arising from the coupled boson-fermion terms. This is 
a key assumption and certainly not traditional as most 
treatments of the LHB focus on the spin physics. How- 
ever, as our emphasis here is on isolating the source of 
zeros in the LHB, demonstrating that the action pos- 
sesses such modes in the absence of the spin-spin scatter- 
ing term would suffice to show that such an interaction is 
indeed ancillary to the essential charge physics. As will 
become evident, our treatment does in fact show this to 
be the case. Under these considerations, the effective low 
energy action turns into 

Sj^ = dr^c^^ [gp{dr - - gte\^] Cka 

- Cq_ktCki + c.c.) I , (2.13) 

where /i denotes the chemical potential, k and q are the 
momenta, and ek = tgije^^'^^^~^^\ This action has a 
BCS-like coupling and hence will have a Green function 
of the form of Eq. (jl.ip . That the Green function must 
be of the form of Eq. (jl.ip is not dependent on the as- 
sumptions delineated earlier. It relies solely on the fact 
that the spectral weight in the lower band exceeds 1 + x 
and hence a new charge e state distinct from the pro- 
jected electrons must be present. Such an excitation can 
only be a composite. 

Consequently, for a given amplitude of (/?q, the 
Fourier transformation of the two point correlator Q = 
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— (7^'^/;^(r)?/^J(0)) becomes 



lUJri 



^ + -(•••), 

M-^tek-i;±(^cJn,k) U 



lUOn 



9t 



(2.14a) 
(2.14b) 

gt/gp, and 



where ujn = (2n + l)7r//3 for n G 
5k, q = 1 — (ck + Cq-k)/^- The ± subscript on the self 
energy arises from the two choices which are possible 
to treat the dynamics of the charge 2e boson. If (fi is 
treated as an independent degree of freedom that can 
condense, then it can be absorbed as a redefinition of 
the interaction strength, s scp. This will correspond 
to a simple condensation of in a non-zero momentum 
particle-particle channel, hence the + sign in front of the 
^t^q-k factor in the denominator of the self energy. As 
will be clear, this is not the interpretation of (p that is 
ultimately consistent with the theory outlined here. Al- 
ternatively, (p^Cia- could be viewed as a new composite 
charge e excitation that results from dynamical spectral 
weight transfer. With such bound modes, the interaction 
term, (/9*c^C|, can be interpreted as a particle- hole scat- 
tering process. To implement this interpretation in the 
Green function, we note that since (p^^c^c^ now describes 
the scattering of an electron Ca- off a composite parti- 
cle (p*VaCa/\^\^ the denominator in the one- loop self en- 
ergy will resemble that of a particle- hole scattering event, 
thereby leading to a — sign in front of the ^tCq-k term 
in the denominator of the self energy. In additon, the 
ellipse in Eq. (|2.14ap represents the terms that originate 
from the mixing between the composite excitations and 
the projected electron, which are at least suppressed by 
the factor t/U. The number gt results from the rescal- 

ing gp^ Cig . Sinc e gt = gt/gp ^ 2x', the t/U 

corrections in Eq. (j2.14ap are comparable to the leading 
term only for < t/2U ^ 0.05. For example, at half- 
filling, the Green function only has the t/U term and the 
spectral weight is governed entirely by the mixing be- 
tween the projected and composite excitations as shown 
previously^^. In the current treatment, we will explore 
entirely the contribution from the leading term which is 
of the form of Eq. (fTT]) . 



D. Free-energy Minimum approach to 2e boson 

Evaluating the Green function is equivalent to a 
random-matrix problem. In the most general case, the 
field cpi must be integrated over with a separate value 
on each site. However, such a mult i- variable integration 
is not tractable in any dimension. From the observa- 
tion that the collective boson is not canonical, that is, it 
does not have its own kinetics, it was previously conjec- 
tured that the spatially homogeneous configuration was 
the most prominent candidate for the ground stat o^^i^^ . 
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FIG. 1. (Color online) Free energy minimization and prob- 
ability distribution of charge 2e boson, (a) Plotted here is 
the AF = F{(poe^^'^) — F{(po) as a function of momentum q. 
(po is the magnitude of the charge 2e boson. As is evident, 
the free energy is strongly dependent on momentum. As ipo 
increases, the minimum of the free energy shifts to (7r,7r). 
This corresponds to a sign change for the bosonic excitation 
around a plaquette. (b) Shown here is the probability distri- 
bution for (po evaluated from the IR action by the relation, 

P{\(p\) = l/Z]^ J D[{c,c}]e~^'^ . The maximum occurs at 
ipQ ^ 1.2 where the corresponding momentum that minimizes 
the free energy is q = (tt, tt). 



Even though such an approach was successful in captur- 
ing some experimental^^ findings, it still leaves an open 
question whether the homogeneous solution minimizes 
the free energy. To this end, we explore some inhomoge- 
neous solutions for (pi to see where the free-energy is a 
minimum. In particular, we explore a staggered configu- 
ration, (pi = e*^"^^|(^o|- It should be noted that a particu- 
lar choice of the configuration of (p does not correspond to 
spontaneous symmetry breaking, since the bosonic mode 
is in lack of inherent dynamics. 

In Fig. [U we directly compute the free energy dif- 
ference between a configuration with a spatial texture 
and the homogeneous state, AF = F{(poe'^^'^) — F{(po), 
where q, determines the spatial dependence of cpi. Ex- 
cept for small values of cpo in which the homogeneous 
solution minimizes the free energy, the distinct minimum 
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occurs at (tt, tt) when the magnitude of the bosonic field, 
(fo^ approaches unity. This is significant because the 
probabihty distribution of (po computed from P{\(p\) = 
J D[{c,c''}]e-^h , has a distinct maximum pre- 
cisely at the value of cp where the (7r,7r) solution min- 
imizes the free energy. This state of affairs obtains be- 
cause a quick inspection of the action reveals that for a 
staggered configuration of (p, the (p^b term actually van- 
ishes. This results in a lowering of the energy relative to 
the homogeneous solution. 

That the (tt, tt) configuration of cpi minimizes the free 
energy is highly significant because the evaluation any 
integral over cpi will be dominated by the staggered so- 
lution. What about the single-particle Green function? 
In our previous treatment of this problem in which we 
assumed that the mixing with the UHB was mediated 
by a homogeneous boson, cpo for all sites, we obtained 
a completely gapped structure at the chemical potential 
for the spectral function. Given that q = (7r,7r) is the 
blobal minimum, our expression for the Green function 
simplifies to 

G{iuJn,k) = J d\(p\\(p\P{(p)g{iu;n,k)\^^^s^^^\^\' 

(2.15) 

The probability distribution function P{(p) is shown in 
Fig. lb. For completeness, we present in Fig. [2] the band 
dispersion corresponding to the maximum in the spec- 
tral function obtained from Eq. (|2.14ap for three different 
cases: 1) homogeneous solution, 2) staggered (tt, tt) phase 
of (pi in the particle-particle channel, and 3) stag- 
gered solution in particle-hole channel, E_. For the ho- 
mogeneous phase, we find a hard gap. Fig. [2fa), because 
no momentum states cross the chemical potential. How- 
ever, as shown earlier, the homogeneous solution does 
not correspond to a minimum in the free energy. Con- 
sider the staggered solutions shown in panels Fig. Efb) 
and Fig. Efc). Fig. Efb) shows that even a staggered 
solution in the particle-particle channel, a gap does not 
occur at the (tt, 0) region of the Brillouin zone. There 
is also a crossing along the zone diagonal. This indi- 
cates that a simple condensation of in a non-zero mo- 
mentum particle-particle channel cannot give rise to the 
no dal/ ant i- nodal dichotomy. There is in fact a clear rea- 
son why cp cannot be treated as an independent degree of 
freedom that can condense. There is a one-to-one corre- 
spondence between Eq. (j2.1Qp and its analogue (Eq. (19) 
of Ref.^^) in the standard perturbative treatment of the 
Hubbard model. In essence, the charge 2e boson replaces 
a string of operators that account for the mixing of double 
occupancy into the lower band. This is why this approach 
is simpler. As it would be completely incorrect to re- 
place that string of operators with an average value, it is 
equally wrong to treat as a variable that can condense. 
In fact, it is well knowr^ that such mean- field trunca- 
tions fail to describe dynamical spectral weight transfer 
in the Hubbard model. 

Consider the third dispersion. Fig. [2fc) in which (p is 
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FIG. 2. (Color online) The low energy band dispersion along 
high symmetry directions for (a) a homogeneous configura- 
tion, (fi = (fo, (b) a staggered configuration (cpi = \(po\e^'^'^'') 
evaluated by Eq. ()2.14ap . and (c) a staggered one evaluated by 
Eq. ()2.14ap . Here, we take the broadening factor 77 = 0.025t 
for a typical value U/t = 10, t^t = -0.3, t'' /t = 0.1, and 
the bare hole doping level as x = 0.12. In the evaluation, the 
parameters, a — x' — x are taken from the numerical esti- 
mates of the number of double occupancies from Liebsch and 



the mediator of a composite charge e state. This corre- 
sponds to a self-energy given by E_. The break in the 
dispersion just above the chemical potential is not fol- 
lowed by a re-entrant crossing at a higher momentum. 
Such a re-entrant crossing would give rise to a closed 
Fermi surface. It is the presence of the additional propa- 
gating degree of freedom which thwarts this re-entrance. 
In addition, there is no crossing at (7r,0), but a broad 
incoherent feature indicative of the pseudogap. Since the 
break-up of the bound state results in a band crossing 
near the (tt, 0) region, the root cause of the pseudogap is 
the bound state formed between the bosonic field, cp and 
a hole as we have advocated previously^. Consequently, 
the pseudogap problem is one of confinement. The corre- 
sponding Fermi surfaces are shown in Fig.[3l The arc-like 
structure is evident. The line of zeros is given by the di- 
vergence of the self-energy and hence it is doublon-holon 
binding that is responsible for killing the intensity on the 
back side of the arc. The Fermi surfaces evolve smoothly 
for the doping levels shown from x = 0.05 to x = 0.18. 
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High 

I 



I 

Low 

FIG. 3. (Color online) The spectral function of the low energy 
theory for each hole doping level, (a) x = 0.05, (b) x = 0.08, 
(c) X — 0.12, and (d) x — 0.18. In the present figures, an 
incoherent background is removed and the intensities of the 
spectral function are normalized for the first quadrant of the 
full Brillouin zone. 

Note also the broad feature at the zone boundary. While 
it is tempting to interpret the broad peak near the antin- 
odal region as an electron pocket, the lack of coherent 
excitations makes this view untenable. 



III. FINAL REMARKS 

The key point this work demonstrates is that two types 
of charge carriers go into forming Fermi arcs. The pro- 
jected electrons are present in any low-energy reduction 
of the Hubbard model and create the spectral weight on 
the high-intensity side of the arc. The zeros correspond 
to composite excitations which are present as a result of 
dynamical spectral weight transfer and hence are present 
only if the UHB is retained or treated appropriately. 
Such composite excitations enter the self-energy through 
the particle-hole channel, as the relevant scattering pro- 



cess is that of a fermion from the composite excitation. 
Both of these features leading to an effective two-fluid 
model^i^ are present within a Wilsonian reduction of 
the high-energy scale in the Hubbard model. The treat- 
ment we have derived here should be valid as long as the 
UHB provides a relevant perturbation to the physics of 
the LHB. Hence, it cannot describe the crossover to the 
Fermi liquid regime in which ip is unbound. Experimen- 
tally, a decoupling of the UHB from the LHB appears to 
take place around x ^ 0.25^^. Accompanying the col- 
lapse is a transition from a small Fermi surface scaling 
with X to a large one with effective area 1—x. The precise 
nature of this transition will be the subject of a future 
study. However, a prediction of this work is that the in 
the pseudogap regime, the volume of the Fermi arc re- 
gion should be given by x' rather than x. This follows 
from the fact that the number of particle-like excitations 
that minimally couple to the electromagnetic gauge is less 
than 1 — x' < 1 — x. Hence, the hole Fermi surface should 
be given by x' . High precision ARPES measurements can 
be employed to verify this result. 

Since our scheme of two type of charge carriers, one 
giving rise to zeros and the other to poles, seems quite 
general, it is tempting to rewrite the IR theory in terms 
of the composite and projected excitations. This would 
require integrating in an additional field for the compos- 
ite, fio- degree of freedom. The composite fermion is not 
canonical, however, and treating it as such would destroy 
the key feature leading to a suppression of the spectral 
on the back-side of the arc. Thus far, we have found no 
consistent way of doing this. Hence, an open problem re- 
mains precisely how the new composite excitation should 
be treated. But that it is present in any standard model 
of Fermi arcs is not in doubt. 
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